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^~~* . Abstract 

OO . 

^5 ' We consider Gaussian ensembles oi m N x N complex matrices. 

We identify an enhanced symmetry in the system and the resultant 
closed subsector, which is naturally associated with the radial sector 
of the theory. The density of radial eigenvalues is obtained in the large 
N limit. It is of the Wigner form only for m = 1. For m > 2, the new 
form of the density is obtained. 
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1 Introduction 

The Wigner semicircle distribution [1] 



/(^) = ~~fy?^^^ ~ ^^' - R < X < R; f{x) =0, \x\ > R , 

TTll 

describes the density of eigenvalues of a gaussian ensemble of single large 
hermitean, symmetric or quaternionic matrices |2j, and finds applications in 
the descripiton of systems in many areas of Physics, from Nuclear Physics to 
Condensed Matter Physics. 

The Gaussian ensemble of a single complex matrix, or equivalently, of 
two hermitean matrices, is also described by a Wigner type distribution of 
eigenvalues of a "radial" variable [3] [1] , the definition of which will be made 
precise in the following. 

It is of great interest and importance to investigate if this continues to be 
a general property of gaussian ensembles of more matrices, and in general, to 
study the properties of systems with a finite number of matrices, particularly 
in their large A^ limit |5]. 

The reason for this importance includes, for instance, the fact that, as it 
has been established already some time ago [6], [7], [8], QCD can be reduced 
to a finite number of matrices with quenched momenta. 

Of more recent interest, the matrix description of D branes [9j has for 
instance lead to the proposal that the large A^ limit of the quantum mechanics 
of the multi matrix description of DO branes provides a definition of M 
theory [10]. In the context of the AdS/CFT duality [H], [IS], [13], due to 
supersymmetry and conformal invariance, correlators of supergravity and 1/2 
BPS states reduce to calculation of free matrix model overlaps [14j, [15] or 
consideration of related matrix hamiltonians p^ . For stringy states, in the 
context of the BMN limit [17] and A/" = 4 SYM, similar considerations apply 
[T8] . [T9] . [20] . A plane- wave matrix theory [21] is related to the M= 4 SYM 
dilatation operator [22] . 

In this communication we will consider m complex N x N matrices 

Za A = l,...,m , 

or equivalently, an even number 2m oi N x N hermitian matrices and a 
gaussian ensemble of such matrices: 



We will obtain the large N description of the system in terms of the 
density of eigenvalues of the matrix 

m 

Y.^aZa (2) 

A=l 

This matrix has a very natural interpretation as a matrix valued radial 
coordinate. 

We will establish that for m = 1 the density of eigenvalues of the radial 
matrix is still described by a Wigner distribution, but that this is no longer 
the case for m >2 and obtain the form of the new eigenvalue distribution. 

In order to obtain these densities, we will first need to establish a new 
result, the measure for the probability density ([I]) in terms of the eigenvalues 
p, = rl t = l,...,Noim- 



A i i>j 

= D^\[dr,rr-'\{rr'-'rr-\r^ -r]f (3) 

i i>j 



m—-l m—-l 
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where the antisymmetric product 

ARM(Pi) = X\p7^ pj^ {pi - Pj) 

i>j 

generalizes the well known Van der Monde determinant A = nj>(Pi ~ Pi)> 
and Cm and Dm are numerical constants. 

This comunication is organized as follows: In Section 2 we restate the 
problem, identify the enlarged symmetries of the gaussiand ensemble and 
identify a complete subset of operators invariant under this enlarged symme- 
try. In Section 3, based on the remarkable fact that this subset of invariant 



operators closes under Schwinger Dyson equations, the Jacobian of the trans- 
formation to these invariant states and to the eigenvalues of ([2]) is obtained. 
In Section 4, the large N density of states on the positive real line is obtained 
for m > 2. In this case, the effective potential contains a new logarithmic 
term which keeps the support of the density of eigenvalues strictly positive, 
i.e., < p_ < p < p+. In Section 5, the density is appropriately extended to 
the whole real line, allowing for a common description of both the m = 1 and 
m > 2 cases. For m = 1, a (restricted) Wigner distribution emerges, whereas 
for m > 2 the required two cut solution is shown to agree with that of the 
previous section, when suitably restricted to the positive real line. In Section 
6, these densities are related to the density of zeros of certain polynomials. 
Section 7 is left for a summary and brief discussion. 

2 Gaussian ensemble and symmetries 

As stated in the Introduction, we consider a gaussian ensemble of m complex 
N X N matrices 

{ZA)ij , A = l,...,m. 
The gaussian potential 

,2 



S, = ^MY.^aZa) 



A 

is invariant under the U{N)"^^^ symmetry 

Za^VaZaV^ A=l,...,m, (4) 

with Va ,A = l,...,m and V unitary matrices. The potential depends only 
on the eigenvalues of the positive definite, hermitean matrix 

E ^aZa , (5) 

A 

which are denoted 

Pi = r'^ , i = 1,...,N , Pi > 0. 

In the context of a physical framework where the m complex matrices 
are naturally associated with 2m hermitean matrix valued coordinates, as 



appropriate, for instance, in the description of the dynamics of branes, these 
eigenvalues have a natural interpretation as the eigenvalues of a matrix valued 
radial coordinate. 

The purpose of this communication is to obtain the large N distribution 
of these eigenvalues for the gaussian partition function 

Z= l\[\[dZA\,dZM,e-''^ 

•^ A ij 

We will do so by "integrating out" the "angular matrix valued degrees 
of freedom or equivalently, by obtaining the jacobian J{pi) of the change of 
variables to the "radial" eigenvalues: 



Z= fl[dp,J{p,)e-'^^'^ 



In order to do so, we will consider correlators of operators that are in- 
variant under the symmetry (jl]), i.e., in the subsector of the theory with this 
enlarged symmetry. Such invariant operators can be built as the trace (sin- 
gle traces) of powers of the matrix (j5]), and hence they depend again on the 
eigenvalues of this matrix only. 

A generating function for such opearators is given by 

$fc = Tre^'^^^^^^^ = ^gifcft = ^e^''^? , 

i i 

or its fourier transform, the density of eigenvalues: 

Hp) = I ^e-^''^^ = Y.^{P- P.) = E ^^P 



r?l 



It turns that these correlators close in this subsector, as it will be shown 
in the following by use of Schwinger Dyson equations 

3 Jacobian 

Schwinger Dyson equations can be obtained from the identity: 

n n dZA\,dZA.,^r^ f ^^^[*]^~'^l = 0' (6) 
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where F[^] is an arbitrary product of invariant operators. This yields: 



< 



92$, 



diZA)l,diZA), 



■F[$] > + < 






> 



<F[$] 



5$* 



dS„ 



d{ZA% ^{ZA)j^ 



> 



(7) 



We denoted in the above: 



< cm >= jwwdZA^.dZM, cm e-'^ = Jmji<^) cm e 



A ij 



Following [23], we now consider the identity 



[rf$] / dk 



Then 



, d 



5$1 






J($) F[^] e 



dk' < 



+ dk' < 



d 



9$, d^y 



d{ZA% ^{ZA),^ 



F[$] > 



d^k d^y 



d{ZA% d{ZA, 



n 






9$ 



+ < 



d<l>k dF[<l>] 

d{ZA% ^{ZA),^ 



> -< F[$] 



k 



dS. 



d{ZA% d^ZA),^ 



>=0 



(8) 



where in the last two terms we used the chain rule. 

Comparing ((Tj) with (jH]), which are equivalent for arbitrary -F[$], it follows 
that: 



dk' 



d{ZA% d{ZA 



+ / dk' 



d 



d^'. 



31 



5$; 



(9) 



d^y 



d{ZA% ^{ZA),^ 
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92$, 



d{ZA%d{ZA),. 



Fourier transforming, and defining 



iLppI 



Wn 



dk r dk' 



27r./ 277 
dk 



Q~ikp^-ik'p' 






-ikp 



92$, 



2vr d{ZA%d{Z, 



tfie differential equation for the jacobian then takes the form 



dp'n. 



91nJ($) 



Vtppi and Wp have been obtained in |4]: 



dp 



/ dVtppi 



d^p') 



(10) 



;ii) 



%p' = dpdp,[p^{p)5{p-p')] 



Wr. 



-d. 



, , , , rfp'$(p') N(m - 1) 



(12) 
(13) 



P-P P 

As a result, the second term in fITT]) vanishes and the Jacobian satisfies: 

'^d^{p) J p~p' P 

This equation was previously obtained in |1], using collective field theory 
methods IMl. The solution is 



In J = / dp^{p)4 dp'^{p) In \p - p'\ + N{m - 1) / dp^{p) In^ 
In terms of the eigenvalues. 



J = IIp^'IIp" p^" \p^ 



m—1 m— 1 

'i ' Pj 



Pj\ = UpT''UpT''pr'(P'^-p^y 



i>j 



Since, up to a constant, J[d^] ~ jYli^Pi^ ^^ have obtained the result 
that we sought to establish: 



/ n n dZA^^JdZM,e-' = / n dp^JiP^) e-"^''^ = [ H "^P^^^P^ 



_ _ -S{P^) 

A ij 

-S{Pi) 



Cm|l[dp^pT 



llpT~'pr'(p^-p^)' 



.i>j 



for potentials invariant under (jl]). 

A couple of comments are in order. When m = 1, the above result reduces 
to the result first obtained in [1], where an explicit parametrization of the 
2N degrees of freedom of two hermitean matrices was obtained, in terms of 
radial and angular matrix valued coordinates. There is a classic result [25] 
[2] parametrizing a single complex matrix in terms of its complex eigenvalues 
and upper diagonal matrix. This parametrization of degrees of freedom, 
useful for holomorphic projections, is different from the one considered in this 
communication. For m = 1, the existence of a closed hermitean subsector 
has also been identified in j26j . 

Gaussian ensembles of rectangular M x N matrices have also been dis- 
cussed in [27| , [2S] and [21] . They can be related to the ensembles discussed in 
this communication when M = mN. In this context, the approach followed 
is equivalent to using the symmetries of the system to set m — 1 of the N x N 
matrices to zero. This corresponds to a "gauge fixed" treatment, as opposed 
to the gauge invariant approach described in this communication, for which 
the eigenvalues r, have a natural identification as a radial coordinate. 

4 Large N density of the eigenvalues 

Writing 

Z = fUdp.Jip,) e-'^^P'^ = fl[dp,e-'^ff , 

i i 

we have 



w"^ 



Seff = —Jdp^p)p (14) 

- dp<l>{p)-l dp'<l>{p')\n\p-p'\-N{m-l) f dp<l>{p)\np , 
I dp^{p) = N . 

In order to exhibit explicitly the N dependence, we rescale p — )■ Np and 
$-)■$, so that 

Seff = N'[—jdp^p)p (15) 

- / dp(^{p)4 dp^{p) In \p - p\ - (m - 1) / dp^{p) Inp] , 

dp^{x) = 1 
The large A^ — > oo configuration is then determined by the stationary 

p-p' A 2p ^ ' 

We note a major difference between the case of one complex matrix (m = 
1) and the case of more than one complex matrix: for more than one complex 
matrix, an additional logarithmic potential is present, in addition to the 
standard Van der Monde repulsion amongst the eigenvalues. As a result, the 
eigenvalues are "pushed away" from p = 0. 

The m = 1 solution will be described in the next section, where the range 
of the density of eigenvalues is extended to the full real line. In this section, 
we obtain the the density of eigenvalues for m > 2. 

The solution to the integral equation ( TT6l) . generalizing that associated 
with Penner potentials [30], is obtained using standard methods [31], together 
with a careful treatment of the p — )• behaviour, following Tan [32]. One 
introduces the function G{z) in the complex plane 



condition dp^ff^ = 0, or 



where $(p) has support only in the interval [p^,p+], p^ > p^ > 0. For large 
\z\, G{z) ~ -, and as z approaches the support of $(p), 

G{p±ni) = i V T^TT'I' p =-r-^ ^T^vr<l>p. 

V P-P 4 2p 

Therefore, this suggests the single cut anzatz 

w^ (m - 1) w;2 ^- 

One also requires [32] that G{z) has no pole as 2; — > 0. These conditions fix: 

2 4 

p± = — (m + l)±— v^ (17) 

It follows that the density of eigenvalues is given by 



111 1^^^^^^^^^^^—^^^^^^^^^^^— 

^(p) = ^V(p+-p)(p-p-) , P-<P<P+ 

= f!./(^-l)(l-^) (18) 

47r y p p 



np V lo w^ 



no longer of the Wigner form. 



5 Symmetric solutions 

In this section, we extend the domain of definition of the density of eigenval- 
ues, allowing us to provide a unified description of the single complex matrix 
case {m = 1) and that of more than two complex matrices {m >2). 
With p = r"^, r > 0, define 

2r<l>(r^) = 0(r) = 0(-r) 

In this way, for an arbitrary function /(r^) 

+00 /"+00 

drfir'')(l)ir) = 2 / dpfipMp). (19) 

00 Jo 
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Returning to (ITB]) . we remark that (p = r^ 



rfp'$(p') p 2r't/r'<l>(r' 



fr'^) 1 /""^ 11 

M = ^^ ^->(0(--7 + --7) 



P — p' Jo r"^ — r 2r 7o 7- _ ^' ^ _|_ ^/ 

2ry_oo ^-^' 
As a result, f lT^ is equivalently written as 

— ^ = — r - ^ ^ , / rfr0(r) = 2 (21) 

When m = 1, this integral equation has the well know Wigner distribution 
as a solution: 



(r) = — -i/^ — r^, —- — <r< 



271 S w^ w w 

We can now write 



47r V w P w'^ 

The symmetric solution of (!2T]) for m > 2 has been discussed in [32] [33] . 
It is a two cut solution, with generating functional 

2 /'^ 1\ „,,2 



w (m — 1) u; 



The cuts are in the intervals [— r_|_, — r_] and [r_, — r+], with r+ > r_ > 0. 
The asymptotic condition and the absence of a pole at z —t- fix 

2 2 4 

r± = ^j(m + 1) ± -^\/m , 

in perfect agreement with ( TT71) . The density is then [33] : 



w 



2 



■a/ (r+^ — r2)(r2 — r^^) , r_^ < ''"^ < 



27r|r| 
This agrees with ( ITSl) . recalling that 2r$(r^) 



11 



r. 



2 



6 Density of eigenvalues and zeros of polyno- 
mials 

It turns out that it is possible to relate the densities obtained in the previous 
section to the density of zeros of certain Laguerre or Hermite polynomials 
[33J L34j. 

We use the results of Calogero's work [35j, based on the classical results 
of Stieltjies [36]. They show that the zeros of the Laguerre polynomial L'^{x) 
satisfy 

A 1 _ 1 / l + a\ 
In terms of eigenvalues, equation (fT6l) takes the form: 



N 



w {m — 1) 



= ^ 23) 

. ...Pi-Pj 4 2pi 

Comparison of these equations shows that the solutions of (l23l) are the zeros 
oiL^-'ifp) ,m>2. 

It has also been established [SS] [3S] that the zeros of the Hermite poly- 
nomial H2N{r) satisfy 

^ 1 
E = ^^ (24) 

j=-N,j^i ^ 

For simplicity, we have chosen an even polynomial. Comparison of the 
above equation with (I2T1) . considered when m = 1 and is expressed in terms 
of eigenvalues, shows that its solutions are the zeros of H2N{,-j^f) [3^ . 

The usual relationship between Hermite and Laguerre polynomials is ob- 
tained, in the discrete, by noting that, with Xi > 0, and since x^j = —Xj, 
the left hand side of S2M can be rewritten as 



N ^ N ^ N 



E ^— = E^— + E -^ 

9^+2^ ^2_^2 ' (25) 



N 
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and comparing with 



7 Summary and discussion 

In this communication, we considered gaussian ensembles of m complex N x 
A^ matrices and identified a closed subsector that is naturally associated with 
the radial sector of the theory. In the large N limit, the ensemble is described 
in terms of the density of radial evenvalues, and these have been obtained: 



m = 1 


m>2 


hp) = i^^, - 1 


^(p) = £v/(^-i)(i-^) 


0<P<^ 


p- <p< p+ 

P± = ^(m + l)±^v^. 


Extending to the full line with <^{p)dp = (j){r)dr , p = r"^, the densities 


e the form: 




m = 1 


m>2 


w — — w 





A (restricted) Wigner distribution is present only for m = 1. 

The existence of this closed sector is related to an enhanced U{N)"^~^^ 
symmetry, and the measure in this subsector has been obtained, with a result 
that generalizes the well known single hermitean matrix Van der Monde 
determinant. 

There are several further areas of study that arise naturally from the 
study presented here. Perhaps of most interest is an investigation of fur- 
ther systmes, such as Hamiltoniand and/or interacting systems, where this 
symmetry is present or where the subsector with this symmetry may pro- 
vide physically relevant truncations. It is also of interest to estabish if the 
results described in this communication extend smoothly to an odd number 
of matrices. 
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